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Abstract 

The one-dimensional dynamics of a classical ideal 'exotic' fluid with equation 
of state p = p{e) < violating the weak energy condition is discussed. Un- 
der certain assumptions it is shown that the well-known Hwa-Bjorken exact 
solution of one-dimensional relativistic hydrodynamics is confined within the 
future/past light cone. It is also demonstrated that the total energy of such 
a solution is equal to zero and that there are regions within the light cone 
with negative (— ) and positive (+) total energies. For certain equations 
of state there is a continuous energy transfer from the (— )-regions to the 
(+)-regions resulting in indefinite growth of energy in the (+) regions with 
time, which may be interpreted as action of a specific 'Perpetuum Mobile' 
(Perpetuum Motion). It is speculated that if it is possible to construct a 
three-dimensional non-stationary flow of an exotic fluid having a finite neg- 
ative value of energy such a situation would also occur. Such a flow may 
continuously transfer positive energy to gravitational waves, resulting in a 
runaway. It is conjectured that theories plagued by such solutions should be 
discarded as inherently unstable. 
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1. Introduction 

If physical laws do not prohibit the presence of exotic matter violating 
the weak energy condition^ and having some other certain properties many 
exciting possibilities arise. For example, solutions of the Einstein equations 
coupled to an exotic matter include wormholes, time machines (e.g. |l|)) 
and cosmologica 



time (see, e 



models with energy density of the Universe growing with 
) for a review and discussion of cosmological consequences) 
leading to the so-called cosmological Doomsday, see e.g. (js). Since theories 
incorporating an exotic matter may lead to counterintuitive and, possibly, 
physically inconsistent effects it appears to be important to invoke different 
thought experiments, which could clarify self-consistency of such theories. 
Here we discuss such an experiment and explicitly show that in a class of 
models containing an exotic matter of a certain kind there could be ever 
expanding with time separated regions of space having positive and negative 
total energies and that absolute values of the energies in these regions could 
grow indefinitely with time while the energy of the whole physical system 
is conserved. This is based on the property of the exotic matter to have 
negative energy density measured by observers being at rest with respect to 
some Lorentz frame and, accordingly, given by the (tt)— component of the 



^Let us remind that the weak energy condition is said to be violated for a matter having 
the stress-energy tensor if there is such a time-hke future directed vector field that 
the inequality T^^f^f^ < holds in some region of space-time. 
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stress energy tensor, T", provided that there are sufficiently large fluid veloc- 
ities with respect to this frame. We also speculate that in a more advanced 
variant of our model there could be an isolated region of space filled by an 
exotic matter with its total energy indefinitely decreasing with time due to 
processes of interaction with some other conventional physical fields carrying 
positive energy. One of such processes could be emission of gravitational 
waves. If conditions for emission of gravitational waves are always fulfilled 
in the course of evolution positive energy is continually carried away from 
the region, which results in a runaway. In this respect it is appropriate to 
mention that the known results on positiveness of mass in General Relativity 
are not valid for the matter violating the weak energy condition, see e.g. (jj) 
and the energy of an isolated region could evolve from positive to negative 
values. Although such a situation resembles the action of a Perpetuum Mo- 
bile of second kind, where heat transfer from a colder part of an isolated 
system to a hotter part occurs, the notion of 'temperature' looks ambiguous 
in our case and, therefore, because of the lack of notation, we refer to this 
hypothetical effect as a 'Perpetuum Mobile of third kind'. In many appli- 
cations a phenomenological description of the exotic matter assumes that 
the matter dynamics can be described by a hydrodynamical model with the 
stress-energy tensor of an ideal Pascal fluid, which can be fully specifled by 
its equation of state, p = p{e), where p is the pressure and e is the comoving 
energy density, which is assumed to be positive below. In this case, violation 
of the weak energy condition can be reformulated as a requirement that the 
pressure is negative with its modulus exceeding the value of the comoving en- 
ergy density. We deal hereafter only with simplest models of hydrodynamical 
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type and neglect the effects of General Relativity and interactions with other 
physical fields. Therefore, in the model explicitly considered in the text the 
transfer of energy from one region of space to another, in particular, between 
the regions having total energies of different signs, is provided by hydrody- 
namical effects. However, as we have mentioned above, it seems reasonable 
to suppose that an analogous runaway effect may happen in a more realistic 
situation, where e.g. energy is continually carried away from a spacial region 
having negative total energy by gravitational waves^. Taking into account 
that gravitational interaction is universal, the proposed 'Perpetuum Mobile 
of third kind' may represent a difficulty in theories, where it emerges. We 
believe that such theories are inherently unstable and must, therefore, be 
discarded. 

It is important to note, however, that in a general scenario the region 

^Note that one should distinguish between the standard hydrodynamical instabihties 
e.g. the ones operating in a fluid having a negative square of speed of sound and the 
instabihty related to violation of the weak energy condition. While the former could lead 
to a highly irregular non-linear dynamics of the system they cannot themselves result in 
formation of spacial regions having a negative total energy, and, a runaway of the kind 
discussed in this Paper, see also the footnote in the next Section. 

Let us also stress that the runaway effect is different from the well-known instabilities of 
linear modes of stationary hydrodynamical flows having a negative energy difference with 
respect to a stationary configuration like the Chandrasekhar-Friedman-Schutz instability 
or the instability operating in shear flows, see e.g. |3) and references therein. 
Although the linear instabilities can lead to a descrease of energy of the main flow they 
cannot bring a system non-violating the weak energy condition to a state with a negative 
total energy. Therefore, in the latter case the runaway effect of the type we consider in 
this Paper is not possible. 
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emitting gravitational waves could have a non-linear three-dimensional dy- 
namics. Explicit solutions of this kind may be quite difficult to obtain due 
to severe technical problems. 



2. Basic definitions and equations 

Let us discuss a one-dimensional planar relativistic flow of an ideal fluid 
with a baratropic equation of state 



where p is the pressure and e is the comoving energy density. As has been 
mentioned in Introduction, we are going to consider later in the text the 
case of an exotic fluid, where the pressure is negative and the weak energy 
condition is violated. For a baratropic fluid this leads to: 



where we introduce the negative of pressure p, a = —p. Since only one- 
dimcnsional flows will be considered, we can also apply our analysis to a 
situation, where a fluid has an anisotropic stress tensor. Say, we can assume 
only one of its components to be non-negative and proportional to S{y,z), 
where y and z are the Minkowski spacial coordinates corresponding to di- 
rections perpendicular to the direction of motion. In this case, equations of 
motion will describe dynamics of a straight string consisting of exotic matter. 

Equations of motion may be written in a divergent form reflecting the 
laws of conservation of energy and momentum 



p = p{e), 



(1) 



a = —p > e, 



(2) 



'tt I rptX ri rptx I rpxx 



(3) 



,t ^ ^,x ^ U> ^^,x = U> 
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where (t, x) are the standard Minkowski coordinates, comma stands for dif- 
ferentiation, and T*-' are the corresponding components of the stress-energy 
tensor: 

T'' = ^\e + v^p), T*^ = 72t;(e + p), T'''' = + vh) , (4) 
where v is the three- velocity and 7 = :^7j=f • 

3. The Hwa-Bjorken solution and the Milne coordinates 



As has been first shown by Hwa (jsl) and later by Bjorken 



the set of 



equations ([3]) has an especially simple 'acceleration-free' solution valid for a 
fluid having an arbitrary equation of state. In this solution velocity of any 
fluid element conserves along the path of the fluid element and the velocity 
field has a very simple form 

v = x/t = e (5) 

For a baratropic fluid the distribution of energy is given by another simple 
implicit relation 

r = exp{- [ , } = exp{- [ y}, (6) 

where 

r = Vt2 _ 3,2^ (7) 

and e=K is a constant of integration. 

Obviously, equations (EHTj) are defined only inside the future/past light 
cone, |r| > |a;|, in an effective two-dimensional Minkowski space described 
by the metric 

ds^ = dt^ - dx^. (8) 



The analytic continuation of the solution on the right/left Rindler wedge 
|r| < |x| is straightforward. 

Although the two-dimensional Minkowski space appears naturally due to 
one-dimensional character of the problem let us remind that the problem 
is defined in four-dimensional Minkowski space. Therefore, for the problem 
with Pascal pressure, where it is assumed that all variables do not depend on 
the coordinates {y, z) perpendicular to x, it is better to say that from the four- 
dimensional point of view the condition |r| = |x| determines four- dimensional 
"light wedges" since it does not depend on directions perpendicular to the 
direction of motion. 

The energy density e is equal to zero on the light cone |r| = if and 
only if the integrals in the exponents in equation ([H]) are positive and diverge 
when e — > 0. Accordingly, in this case, the condition ([2]) must be satisfied. 
Additionally, in order to make the integrals divergent we must have 



when e 0. Provided that condition ([9]) is valid the solution may be con- 
sidered as confined within the future/past light cone. For simplicity we are 
going to consider only this case in our analytical calculations later on^ . 
The solution ([MS]) has a self-evident form in the Milne coordinates (r, y), 

■^In this case, the corresponding hydrodynamical models are unstable with respect to 
growth of small perturbations. One may proceed, however, either assuming that these 
models are valid only for the considered types of hydrodynamical flows or considering 
them as effective models invalid for sufficiently large perturbation wavenumbers. In any 
case, we expect that our main conclusions do not depend on whether the considered models 
are hydrodynamically unstable or not. 




(9) 
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where the time r is defined by equation ([7j) and we introduce the rapidity 



?/ = lnWl^ (10) 

as a new spacial coordinate. In these coordinates the metric (IHl) has the form 

ds^ = dr^ - T^dy\ (11) 

Although the metric (fTTi) may be understood as describing an expanding 
one-dimensional spacially-uniform universe with the scale factor a(r) = r, 
obviously it corresponds to the same flat Minkowski space since it is obtained 
from the metric (IHl) by the coordinate transformation (JTj), (|TOl) . see e.g. Iiol ) 
for an additional discussion of the Milne coordinates. Taking into account 
that the transformation law between the velocity v deflned with respect to 
Minkowski coordinates (t, x) and the peculiar velocity v deflned with respect 
to the orthonormal frame associated with the Milne coordinates is determined 
by the relation 

1 + 4f 

it is clear that the solution ([Mn]) is simply the spacially-uniform solution in the 
Milne coordinates with the peculiar velocity -i; = . In particular, equation 
([6]) immediately follows from the flrst law of thermodynamics written for 
an adiabatic expansion of a fluid having distribution of its thermodynamical 
variables uniform with respect to the coordinate y: 

where \^ oc r is a comoving volume. 
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4. Properties of the solution 

4.I. The energy integral 

Provided that the condition Q is valid we can assume that the energy 
density and, accordingly, the components of the stress-energy tensor are dif- 
ferent from zero only within the future/past light cone. For simplicity, let us 
consider only the future light cone implying that t > from now on. In this 
case, equations ([3]) yield that the total energy of the flow 

E = j dxT^* = t j d^T^* (14) 

does not depend on time t. 

Let us show that this integral is precisely equal to zero for an exotic 
baratropic fluid satisfying the condition (Q. Taking into account that the 
distribution of energy density and pressure are even functions and the velocity 
distribution is an odd function of the spacial coordinate x, respectively, it 
suflices to prove that the quantity 

is equal to zero. From equation ([6]) it follows that 

de , , 

a = T— + e 16 
dr 

and, therefore, we have 

Now we change the integration variable from ^ to r keeping the value of 
t flxed. Taking into account that 
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T 



(19) 



we obtain 
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1 r erdr 



(20) 



tj VW^' 



where the values of r and e corresponding to ^ = 0, 1 are omitted. Integrating 
by parts the first integral in ( 120|) and taking into account that the boundary 
terms are equal to zero for the solution ([MSD satisfying the condition ([9]), we 
obtain 



4-2. Lorentz invariance and vacuum-like nature 

It is easy to see that solution (IMHl) has the same form in all coordinate 
systems connected by the Lorentz transformations: (t,x) {t',x'). Indeed, 
as follows from equation ([7j) the time r is invariant under the Lorentz trans- 
formations. Therefore, equation ([6]) contains only invariant quantities and 
is the same in all Lorentz frames. It is also evident that when the Lorentz 
transformations are considered the quantity ^ = x/t and the three- velocity 
V are transformed in the same way. Therefore, from equation ([5]) it follows 
that the same equation is valid for the transformed quantities. 

It is clear that the total energy and momentum of the flow are equal 
to zero in all Lorentz frames. This is frequently considered as being the 
definition of vacuum solutions in different theoretical schemes. Thus, one 
may state that solution (JME]) plays a role of a non-trivial vacuum solution 
for the exotic fluids satisfying ([9]). 




(21) 
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4-3. A hypothetical model of Perpetuum Mobile 

As has been mentioned in Introduction, the very possibihty of existence 
of solutions having negative/zero total energy is determined by the fact that 
for equations of state violating the weak energy condition the energy density 
determined with respect to a fixed Lorentz frame 

oc e - v^a (22) 

can be negative provided that the fiuid velocity is sufficiently large, 

V > Vcrit = \ -■ (23) 



In the case of our solution the space bounded by the light cone condition |x| < 
t is divided into a set of regions having opposite signs of T** and, accordingly, 
different signs of the total energy. In what follows let us refer to the regions 
with T** > (T** < 0) as (+)-regions ((— )-regions). The coordinates of 
boundaries between the (+) and (— )-regions, Xcrit, can be found from the 
condition v = Vcrit and, respectively, from the implicit equation 



where Tent = ^/t^ — xl^u- In general, equation (!24|) could have several roots 
on the interval < x < t with values depending on equation of state. How- 
ever, for a reasonable equation of state with sufficiently smooth dependence 
of 0" on e there must be a (— )-region adjacent to the light cone boundary 
X = t and a (+)-region close to the point a; = 0. The total energy of the 
region adjacent to the light cone, E_, can be easily calculated from equa- 
tions ffn|) . ffT^ . and (1201) . where we take into account that after integration 
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by parts of equation (12UI) only the boundary term at x = Xcrit contributes to 
the integrah 

E-= f T''dx = -Xc„te, (25) 

where Xcrit denotes the largest root of equation fl2^ in the interval < x < t 
from now on. Taking into account the fact that the total energy is zero and 
using the symmetry of the problem, we see that the energy in the interval 
< x < Xcrit, -E+ = —E^, and therefore 

= / r'dx = ^ = exp ( / (26) 

where we use equations ([HD, (1^ . (12^ and all quantities are assumed to be 
evaluated along the world line determined by equation ^IM . It is instructive 
to calculate the time derivative of -E+ differentiating the integral in (l26l) on 
time and using equations ([3]) and (1231) to obtain: 

E^ = av = ^fae. (27) 

From equation (1271) it follows that the energy of the region < x < Xcrit 
constantly grows with time due to energy flow from the (-)-region Xcrit < x < 
t. In principal, there could be two possibilities for the asymptotic behaviour 
of -E+ in the limit t ^ oo depending on form of the equation of state: 1) there 
could be a finite asymptotic value of E^, and 2) a finite asymptotic value 
could be absent and the energy Ej^ could infinitely grow with time. The latter 
case represents a specific instability, where there is infinite growth of energy in 
one region of space and infinite decrease of energy in the other region. Let us 
suppose that physical laws do not prohibit existence of the hydrodynamical 
systems violating the weak energy condition and having solutions of this type. 
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Assuming that such a system could be made by an advanced civihzation, 
which could also ensure that utilization of energy released in the region < 
X < Xcrit does not significantly perturb the solution this could provide an 
infinite source of energy. Therefore, this type of solution may be classified as 
a hypothetical 'Perpetuum Mobile'. Such a solution is discussed below. 

4.4- An explicit example 

Let us specify the equation of state and consider the simplest possible 
case of linear dependence of cr on e : 

cr = (l + a)e, (28) 

where the parameter a > 0.^ In this case integration of equation ([6]) gives 

e = Cr", (29) 

where the parameter C > 0. The critical velocity Vcrit = "y^j a.nd, accord- 
ingly , 

Xcrit = (30) 

V 1 + a 

From equations d?]), ([24]), ([27D, ([28D and ([29D we obtain 



^Let us stress that there is no well defined limit to the case of the 'standard' vacuum 
equation of state a = (i.e. p — — e) and that all our conclusions are not applicable to 
this case. 

^Note that in this case it is obvious that the hypersurfaces separating the (+) and (— ) 



regions are time-like. 
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From equation fl3ip it follows that in the case of the linear equation of state 
there is one (+)-region and one (— )-region in the range < x < t. The energy 
of the (+)-region ((— )-region) grows (decreases) indefinitely. Therefore, the 
simplest linear equation of state determines solution, which may be classified 
as the 'Perpetuum Mobile'. 

5. Discussion 

When the weak energy condition is violated the Hwa-Bjorken solution is 
likely to be not the unique one having the total energy E equal to zero. Say, 
in the model with the linear equation of state fl28p it is easy to find a family 
of self-similar solutions, where the velocity f is a function of the self-similar 
variable f = /i(^), and the energy density e has the form: 

e = t^f2{0- (32) 

Substituting these expressions in equations ([3]) and assuming that the equa- 
tion of state is given by (125]) one can easily get two ordinary differential 
equations for the functions /i and f2- In this case we have 

/ + 00 P + OO 

dxT'' = t^^' «(0, (33) 
-oo J ~oo 

where /s is expressed through /i and f2- It is assumed that the energy 
density is either equal to zero when ^ is sufficiently larger or tends to zero 
with increase of ^ fast enough to make the integral convergent. When this 
condition is fulfilled and (3^—1 the energy E is equal to zero. Indeed, the 
energy must not depend on time. On the other hand it is seen from fl33l) that 
the energy is proportional to t'^^^. This means that the integral dC,/^^^) 
must be equal to zero provided that (3 ^ —1. 
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A more difficult and interesting problem would be to construct an explicit 
solution having a negative value of the total energy. In this case one should 
invoke more sophisticated methods of one-dimensional relativistic hydrody- 



namics such as e.g . the hodograph method introduced by Khalatnikov (I 111 ) 
(see also (Il2l ). (Il3l)). In the case of the linear equation of state fl2Sl) it would 
also be interesting to exploit the formalism developed in Ref . ( jl4| ) , where the 
set of equations ([3]) is reduced to a single equation, which can be analysed 
for new solutions. 

Another approach to the problem consists in use of numerical methods. 
It seems that numerical methods are more suitable for hydrodynamically 
stable models e.g. based on the 'Chaplygin type' equation of state p = — e^/e, 
where e^, is a constant. In framework of numerical methods one can consider 
a hydrodynamical motion with fixed boundary conditions on a fixed spacial 
interval, e.g. v = when x = and when x = Xi, where Xi is fixed. The 
total energy of the motion determined by initial conditions can, in principal, 
be either negative or zero for fluids violating the weak energy condition. 

An interesting development of these studies would be to consider a model, 
where the exotic matter having anisotropic tension is concentrated on a 
straight one-dimensional line in three dimensional space. One-dimensional 
motion excited on the line could produce gravitational waves carrying away 
positive energy from the system. Accordingly, the energy of the motion could 
decrease indefinitely provided that generation of gravitational waves persists 
in the course of evolution of the system. Being extrapolated on realistic 
three-dimensional motions this effect could be dangerous for models of exotic 
matter, which can be effectively described in the hydrodynamical approxi- 
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mation. The energy of such motions could decrease indefinitely, resulting in 
a runaway. In this case, the corresponding models should be discarded. 

It would be also interesting to look for a quantum model having the 
properties disscussed in the Paper. Note that in quantum case some exotic 
properties of behaviour of field systems may be expected even when they 



have a "normal" classical limit. Say, as was discussed e.g. in Ref (jlSl), in 



the Milne universe vacuum expectation value of the comoving energy density 
can be negative even for the simplest model of a non-interacting scalar field 
with sufficiently small mass, in a certain vacuum state. 

At the end we would like to point out that although the runaway pro- 
cess related to emission of gravitational waves may be technically difficult to 
construct for the exotic matter with positive comoving energy density, it can 
be easily constructed for even more exotic 'ghost' matter having a negative 
value of the energy density in all frames implemented by physical bodies and 
clocks. For example, we can make use of the model of rotating relativistic 
string with two monopoles at its ends emitting weak gravitational waves, see 
Ref. m. It is sufficient to change the sign of the Lagrangians describing the 
string and the monopoles while keeping the sign of the gravitation part of the 
action fixed to convert this model to a model of a 'ghost' matter interacting 
with gravity. It is clear that neither the string dynamics nor conditions of 
emission of gravitational waves are significantly affected by this procedure. 
Such a model describes a finite length string with its length ever increasing 
with time thus making the total energy of the string-monopoles system ever 
decreasing. The positive energy carried away by gravitational waves may be 
exploited by an advanced civilization able to constract such a device. 



16 



I am grateful to V. I. Dokuchaev, I. D. Novikov, V. N. Lukash, G. W. 
Gibbons, J. C. B. Papaloizou and Vovik Strokov for useful remarks. It is my 

sincere pleasure to thank D. A. Kompanccts for many important discussions 
and comments. This work has been supported in part by RFBR grants 
07 - 02 - 00886-a and 08 - 02 - 00159-a and by the governmental grant 
NSh-2469.2008.2. 

References 

[1] M. S. Morris, K. S. Thorne, U. Yurtsever, Phys. Rev. Lett., 61, 1446 
(1988) 

[2] E. J. Copeland, M. Sami, S. Tsujikawa, Int. J. of Mod. Phys. D, 15, 
1753 (2006) 

[3] R. R. Caldwell, M. Kamionkowski, N. N. Weinberg, Phys. Rev. Lett, 
91, 071301 (2003) 

[4] R. Schoen, S. T. Yau, Phys. Rev. Lett. 43, 1457 (1979); R. Schoen, S. 
T. Yau, ibid 48, 369 (1982); G. T. Horowitz, M. J. Perry, ibid 48, 371 
(1982) 

[5] S. Chandrasekhar, Phys. Rev. Lett., 24, 611 (1970) 
[6] J. L. Friedman, B. F. Schutz, Ap. J., 222, 281 (1978) 
[7] W. Glatzel, Reviews in Modern Astronomy 4, 104 (1991) 
[8] R. C. Hwa, Phys. Rev. D 10, 2260 (1974) 
[9] J. D. Bjorken, Phys. Rev. D 27, 140 (1983) 

17 



[10] N. D. Birrell, P. C. W. Davies, "Quantum Field in Curved Space", 
Cambridge University Press, Cambridge, UK (1982) 

[11] I. M. Khalatnikov, Zh. Eksp. Teor. Fiz. 27, 529 (1954) 

[12] E. P. T. Liang, Astrophys. J. 211, 361 (1977) 

[13] P. Carbonaro, Phys. Rev. E 56, 2896 (1997) 

[14] A. Bialas, R. A. Janik, R. Peschanski, Phys. Rev. C 76, 054901 (2007) 

[15] T. S. Bunch, S. M. Christensen, S. A. Fulling, Phys. Rev. D 18, 4435 
(1978) 

[16] X. Martin, A. Vilenkin, Phys. Rev. D 55, 6054 (1997) 



18 



